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Abstract 

Using the eigenvalues and eigenvectors of correlations matrices of some of the main financial market indices 
in the world, we show that high volatility of markets is directly linked with strong correlations between them. 
This means that markets tend to behave as one during great crashes. In order to do so, we investigate 
financial market crises that occurred in the years 1987 (Black Monday), 1998 (Russian crisis), 2001 (Burst 
■ of the dot-com bubble and September 11), and 2008 (Subprime Mortgage Crisis), which mark some of the 

| largest downturns of financial markets in the last three decades. 
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1 — '; 1 Introduction 

<n : 

J> , The study of why many world financial markets crash simultaneously is of central importance, particularly 
0^ ' after the recent worldwide downturn of the major markets in 2007 and 2008. Economists have been studying 

■ the reasons why markets crash, and why there is propagation of volatility from one market to another, since a 
t— i . long time. After the crash of 1987, many studies have been published on transmision of volatility (contagion) 
£\j | between markets using econometric models [TJ-[I5], on how the correlation between world markets change 

■ with time [E]-[T9], and how correlation tends to increase in times of high volatility |20J-[34J. This issue is of 
particular importance if one wishes to build portfolios of international assets which can withstand times of crisis 

• • | [35] -[47]. Many models were proposed by both economists and physicists in order to explain the correlation of 
.£h ' international financial markets [15]-[BB], which is considered a complex system with many relations which are 
^ . difficult to identify and quantify. 

One tool that was first developed in nuclear physics for studying complex systems with unknown correlation 
structure is random matrix theory [67J-[70j, which confronts the results obtained for the eigenvalues of the 
correlation matrix of a real system with those of the correlation matrix obtained from a pure random matrix. 
This approach was successfully applied to a large number of financial markets |71j-[100j, and also to the relation 
between world markets [101j-[102j. This approach was also used in the construction of hierarchical structures 
between different assets of financial markets [103j-[142j. 

Our work uses the tools of random matrix theory to analyze the correlation of world financial markets in 
times of crisis. In order to do so, we use data from some of the largest worldwide crashes since 1980, namely 
the 1987 Black Monday, the 1998 Russian Crisis, the Burst of the dot-com bubble of 2001, the shock after 
September, 11, 2001, and the USA subprime mortgage crisis of 2008. We start by defining a global financial 
crisis (Section 2) based on evidence of some financial markets chosen from diverse parts of the world. Then, 
we discuss some of the main theoretical results on Random Matrix Theory (Section 3). Then, in Section 4, we 
make a quick discussion on how we collected the data and how it was treated. 

In the next three sections we study the correlation matrices between the log-returns of a number of financial 
market indices chosen so as to represent many geographical parts of the world and a diversity of economies 
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(sections 5 to 8). In each section, we calculate the eigenvalues of the correlation matrix of the chosen indices 
and then study the eigenvector that corresponds to its largest eigenvalue, which is usually related with a market 
mode, which is a co-movement of all indices. We then calculate correlation matrices in running windows and 
compare the average correlation between markets with the volatility and average volatility of the market mode 
obtained previously, showing that times of large volatility are strongly linked with strong correlations between 
world financial indices. 

Section 9 looks more closely at the probability distribution of the correlation coefficients in different intervals 
of time and tests the hypothesis that it becomes closer to a normal during periods of crises. 

Since the study of world stock exchanges involve dealing with different operating times, in Section 10, we 
compare the results obtained in the previous sections with results obtained by using the log-returns of Western 
markets with the log-returns of the next day in Asian markets. We also compare the results obtained in the main 
text of the article with those obtained by using Spearman's rank correlation instead of Pearson's correlation. 

We finish by discussing recent methods for studying random matrices obtained from non-Gaussian distri- 
butions, shuch as t-Student distributions, which represent more closely the probability density distributions 
obtained from financial data [145]- [US] , in assossiation with other measures of co-movement of financial indices 
that are more appropriate for systems with very strong correlations, like it happens in times of financial crises 
[117] [H5] . 

Since world stock market indices (countries) are easier to relate with than equities in a stock market 
(companies), one of the aims of this article is to be a pedagogical introduction to most of the techniques that 
are used when Random Matrix Theory is applied to financial data. Hence we also supply an ample bibliography 
on the subject. 



2 Defining a global financial crisis 

Before studying periods of financial crises, we must make it clear what we consider to be a global crash of the 
financial markets. In order to adopt a more precise definition, we considered the time series of 15 financial 
markets representing different regions of the world from the beginning of 1985 until the end of 2010. Looking 
at the closing indices of every day in which there was negotiation, we considered the log-returns, given by 

St = ln(P t ) - ln^) « Pt ~ Pt - 1 , (1) 

what makes it easier to compare the variations of the many indices. After that, the 10 most negative variations 
were chosen. 

In order to illustrate the procedure, we consider the Dow Jones index of the New York Stock Exchange 
(NYSE). Figure 1 shows the log-density distribution for this index with data from 01/02/1985 to 12/51/2008. 
The log-density, defined as 

log-density = ln(l + density) , (2) 
is used instead of simple density in order to better visualize the most extreme points. 
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Figure 1: log-density distribution of the Dow Jones index of the NYSE, from 01/02/1985 to 12/31/2008. 
The ten most negative values of the log-returns are bellow —0.07. These events occured in the following 
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occasions: 10/19/1987 (22.61%), 10/26/1987 (8.04%), 01/08/1988 (6.85%), 10/13/1989 (6.90%), 10/27/1997 
(7.18%), 09/17/2001 (7.13%), 09/29/2008 (6.98%), 10/09/2008 (7.33%), 10/15/2008 (7.87%), and 12/01/2008 
(7.70%). These dates include the 1987 Black Monday, part of the Asian Crisis of 1997, the 1998 Russian Crisis, 
the aftermath of September 11, 2001, and the Subprime Mortgage Crisis of 2008. 

The same technique was used for the Nasdaq (USA), S&P/TSX Composite from Canada (Can), Ibovespa 
(Brazil), FTSE 100 (UK), DAX (Germany), ISEQ (Ireland), AEX (Netherlands), SENSEX 30 (India), Colombo 
All-Share (Sri Lanka), Nikkei (Japan), Hang Seng (Hong Kong), TAIEX (Taiwan), Kospi (South Korea), Kuala 
Lumpur Composite (Malaysia), and Jakarta Composite (Indonesia). The next table displays the years of major 
drops (between the beginning of 1985, and the end of 2010) and the quantity of markets which presented those 
falls. When a market drops substantially more than once in the same year, these are counted more than once 
as well, in order to gauge the depth of the shocks. This helps identity the times where there were major crashes 
around the world. 
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Year 


1998 


1999 


2000 


2001 
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2010 


Occurences 


11 


1 
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1 


3 
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3 








50 


2 






Table 1: number of occurences per year of major drops in fifteen diverse stock markets in the world. 



It is possible to pinpoint two major crises in 1987 and 2008, and minor crises in 1989, 1990, 1997, 1998, and 
in 2001. The crisis of 1987 corresponds to the so called Black Monday, the one in 1989 is the USA saving and 
loan crisis, 1990 are the Japanese asset price bubble and the Scandinavian banking crisis, 1992 is the so-called 
Black Wednesday, 1997 is the Asian financial crisis, 1998 is the Russian crisis, 2000 and 2001 correspond to the 
Burst of the dot-com bubble, and 2008 corresponds to the Subprime Mortgage Crisis in the USA. 

We shall apply a theory called Random Matrix Theory in order to analyze four of these crises. The next 
section gives a pragmatic introduction to this theory. 



3 Random matrix theory 

Random matrix theory had its origins in 1953, in the work of the German physicist Eugene Wigner |67| |68| . He 
was studying the energy levels of complex atomic nuclei, such as uranium, and had no means of calculating the 
distance between those levels. He then assumed that those distances were random, and arranged the random 
number in a matrix which expressed the connections between the many energy levels. Surprisingly, he could 
then be able to make sensible predictions about how the energy levels related to one another. 

This method also found connections with the study of the Riemann zeta function, which is of primordial 
importance to the study of prime numbers, used for coding and decoding information, for example. The 
theory was later developed, with many and surprising results arising. Today, Random Matrix Theory is 
applied to quantum physics, nanotechnology, quantum gravity, the study of the structure of crystals, and may 
have applications in ecology, linguistics, and many other fields where a large amount of apparently unrelated 
information may be understood as being somehow connected (for a recent book on the subject, see |70|). The 
theory was also applied to finance in a series of works dealing with the correlation matrices of stock prices, and 
also to risk management in portfolios [71j-[100 j (for recent reviews on the subject, see |93| and |94|). 

In this section, we shall focus on the results that are most important to the present work, which is studying 
the correlations between world financial markets in times of crisis. The first result of the theory that we shall 
mention is that, given anlxJV matrix with random numbers built on a Gaussian distribution with average 
zero and standard deviation a, then, in the limit L — > oo and N — > oo such that Q = L/N remains finite and 
greater than 1, the eigenvalues A of such a matrix will have the following probability density function, called a 
Marcenku-Pastur distribution |69| : 

Q V(A + -A)(A-A_) 

pw ~ w x ' (3) 
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where 

Q ~ 2 V^) 

and A is restricted to the interval [A_, A+]. 

Since the distribution (jSJ) is only valid for the limit L — > oo and N — > oo, finite distributions will present 
differences from this behavior. In figure 2, we compare the theoretical distribution for Q = 10 and a = 1 to 
distributions of the eigenvalues of three correlation matrices generated from finite L x N matrices such that 
Q = L/M = 10, and the elements of the matrices are random numbers with mean zero and standard deviation 
1. 
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Figure 2: theoretical and sample finite distributions for a random matrix with N = 10 and L = 100, 
N = 30 and L = 300, and TV = 100 and L = 1000. 

Consequently, real data will deviate from the theoretical probability distribution. Nevertheless, the theo- 
retical result may serve as a parameter to the results obtained experimentally. 

Another source of deviations is the fact that financial time series are better described by non-Gaussian 
distributions, such as t Student or Tsallis distribution. This can be seen from figure 1: a Gaussian distribution 
would be represented by a parabola, what is clearly not the case. Recent studies |143| - [T4"8] developed part of 
the theoretical framework in which finite series and series with fat tales, as is the case of financial time series 
of returns, can be studied. 

Since Random Matrix Theory is based on random matrices with a single standard deviation a, we must 
compensate the data obtained from the many indices so that all series have average zero and the same standard 
deviation, which we chose to be equal to one. This can be done using the formula 



s t - (S) 



a 



(5) 



where (S) is the average of the time series used, and a is its standard deviation. 



4 Data 

We shall work with one stock market index of each country (with the expeeption of the USA, with two indices). 
The indices were chosen among the ones that are considered benchmarks in each stock market. The data were 
collected from Bloomberg, from 1980, or the first available date, until the end of 2010. For 1987, we collected 
23 indices (4 from the Americas, 9 from Europe, and 10 from Asia); for 1998, we have 63 indices (12 from 
America, 24 from Europe, 19 from Asia, 1 from Oceania, and 7 from Africa); from 2001, 79 indices (13 from 
America, 29 from Europe, 26 from Asia, 2 from Oceania, and 9 from Africa); for 2008, we have 92 indices (14 
from America, 35 from Europe, 30 from Asia, 2 from Oceania, and 11 from Africa). The number of indices 
collected grew in time both due to the adoption of the indices by their respective stock markets, teh availability 
of data, and by the emergence of new countries. 

This work is motivated by the will to understand how each index affects the others, so as to later attempt 
to build a model of how crises propagate in a network of indices. Thus, we do not consider here the indices 
normalized to a single currency, which would be useful, as an example, for building portfolios for investors. 
That is because we want the numbers to be the ones to which agents operating in their own stock markets 
react when they take decisions. That is also the reason we are not using indices that are standardized in terms 
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of methodology, like the ones calculated by Morgan and Stanley Capital International (MSCI), which are used 
mainly by researchers and international investors, but that are not the ones usually published by the press, or 
seen on the news broadcasts. One of the authors (LSJr) has some research underway using MCSI indices. 

When analyzing the data, we had to be careful with the differences in public holidays or weekends among 
countries. Particular care had to be taken with Israel, Palestine, Jordan, Saudi Arabia, Kwait, Bahrein, Qatar, 
the United Arab Emirates, Ohman, Bangladesh, and Egypt, for which weekends do not occur on Saturdays 
and Sundays, but on Fridays and Saturdays, for example. For these countries, we shifted data so as not to 
lose information, making missing data due to weekends coincide with the other indices. Our general rule was 
that, when more than 30% of the markets didn't open on a certain day, we removed that day from our data, 
and when that number was bellow 30% , we kept the existing indices and repeated the last computed index 
for each of the remaining ones. We did not make linear extrapolations of missing indices, for we could then 
lose the effects of drops, like the one that occured after September 11, when the stock exchanges of the USA 
remained closed for some days. 

Another pressing problem was that markets do not operate at the same time zones, so we had to decide 
wether to consider the data concerning American countries at the same day as Asian-Pacific countries, for 
example, or to shift the data for Asian countries so as to compare indices from the USA with the next day 
index from Japan, for example. There is even some evidence [27] that the correlations of Asian with the USA 
indices increase when one considers the correlation of the USA indices with the next day indices of the Asian 
markets. 

We decided to consider all indices taken at the same date. This was motivated by some comparisons 
between the correlations among indices: we compared the individual correlations among most of the indices 
and checked if the correlation increased or decreased shifting the relevant data. The result was inconclusive, 
for there was an almost equal number of correlations which increased taking Western indices one day before 
their Eastern counterparts as there were other correlations that remained higher taking the same day indices 
(some correlations were higher when we took the Western indices one day after the Eastern ones). In order to 
gauge the effects of such a shift, we recalculated all our results by shifting the indices from Asia and Oceania 
in one day, without many changes in the outcomes. Those results are better explained in Section 10. 

One option, frequently adopted, to avoid the problem of different operating hours between international 
markets is to consider weekly data instead of daily data. We didn't adopt that approach so as not to miss 
major changes in markets, which tend to occur during a small interval of days. Instead, we preffered to compare 
results with and without shifting part of the data by one day. 

In the next four sections, we shall consider the years 1987, 1998, 2001, and 2008 in the light of Random 
Matrix Theory. The indices we used, together with their countries of affiliation, the symbols used for them, 
and their codes in Bloomberg, are placed in table 2, in the Appendix. 

5 1987, the Black Monday 

In 1987, the financial world lived a time of panic, much like the one of the great crash of 1920. In a matter of 3 
days, most of the stock markets in the world lost about 30% of their value and trillions of dollars evaporated, 
leaving a trace of destruction that affected what is reffered to as real economy for many years. The day the 
first and major collapse occured, a Monday, was later called the Black Monday. 

In order to analyze that crisis, we consider now the correlation matrix of 23 indices of stock exchanges 
around the world: the S&P 500 from the New York Stock Exchange (S&P), and NASDAQ (Nasd), both from 
the USA, S&P TSX from Canada (Cana), Ibovespa from Brazil (Braz), FTSE 100 from the United Kingdom 
(UK), ISEQ from Ireland (Irel), DAX from Germany (Germ), or West Germany in 1987, ATX from Austria 
(Autr), AEX from the Netherlands (Neth), OMX from Sweden (Swed), OMX Helsinki from Finland (Finl), 
IBEX 35 from Spain (Spai), ASE General Index from Greece (Gree), SENSEX from India (Indi), Colombo All 
Share from Sri Lanka (SrLa), DSE General Index from Bangladesh (Bang), Nikkei 25 from Japan (Japa), Hang 
Seng from Hong Kong (HoKo), TAIEX from Taiwan (Taiw), Kospi from South Korea (SoKo), Kuala Lumpur 
Composite from Malaysia (Mala), JCI from Indonesia (Indo), and the PSEi from the Philipines (Phil). So, we 
have three indices from North America, one from South America, nine from Europe, and ten from Asia, with a 
total of 23 indices. These offer a good variety of indices worldwide. In subsequent years, we shall increase this 
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number, mainly due to the appearance of new indices and countries, and to the access to data about them. 

We shall use 1987 as an example for the other years, and because of that we will be showing more details 
in the calculations for that year. Calculating the correlation matrix for the indices that are being considered, 
one obtains a 23 x 23 matrix. The average of the values of this correlation matrix is a good measure of the 
overall correlation between the many indices (the average is taken over the elements of the correlation matrix 
for which i < j). For the present correlation matrix, the average is given by < C >= 0.16, with standard 
deviation a = 0.04. Since the correlation matrix is 23 x 23, and the number of days considered in calculating 
it is 256, we then have Q = L/M = 256/23 ~ 11.130, and the upper and lower bounds 



A_ = 0.490 



A j 



1.689 



(6) 



for the eigenvalues that constitute the bulk of the eigenvalue distribution due to noise. 

A frequency distribution of the 23 eigenvalues of the correlation matrix is shown in figure 3, with the 
theoretical distribution of an infinite random matrix for Q = 11.130 with mean zero and standard deviation 
one superimposed on it. In figure 4, the eigenvalues are plotted in order of magnitude. The shaded area 
indicates the region predicted by theory were the data related with a purely random behavior of the normalized 
log-returns. 

P(A) 

1.2 -- 



1 -- 



0.8 



0.6 -- 



0.4 -- 



0.2 -- 




1 1 1 1 — L- 1 -^ A 

1 2 3 4 5 6 

Figure 3: frequency distribution of the eigenvalues of the 
correlation matrix for 1987. The theoretical distribution 
for a random matrix is superimposed on it. 
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Figure 4: eigenvalues in order of magnitude. The 
shaded area corresponds to the eigenvalues predicted 
for a random matrix. 



Only 60% of the eigenvalues fall inside the region predicted by Random Matrix Theory. Note that the 
highest eigenvalue stands out from all the others, being more than three times bigger than the uppermost limit 
A + of the theoretical distribution. This is in agreement with many other results, obtained for a great number 
of financial institutions to which this same formalism has been already applied [71] - [100 ]. It is believed that 
this eigenvalue corresponds to the action of a single market, which influences all the other members of the 
correlation matrix. Usually, for the correlation matrix of individual stocks in a single market, this eigenvalue is 
much larger, some times 25 times larger, than the largest eigenvalue predicted for the correlation matrix of a 
random time series, although the size of the sample directly influenciates that as well. In our case, it responds 
for about 28% of the collective behavior of the time series being considered, which is the ratio of the largest 
eigenvalue and the sum of the eigenvalues of the correlation matrix. 

Figure 5 shows the contributions of the many indices which we are considering in our study in some of 
the eingenvectors of the correlation matrix. The blue bars represent positive values and the red bars represent 
negative ones. 



t 1 






1 



S&P Nasd Cana Braz UK Irel Germ Autr Neth Swed Finl Spai Gree Indi SrLa Bang JapaHoKoTaiw SoKo Mala Indo Phil 



6 



0.6 -- 




S&P Nasd Cana Braz UK Irel Germ Autr Ncth Swed Finl Spai Gree Indi SrLa Bang Japa HoKoTaiw SoKo Mala Indo Phil 



e23 



0.6 + 
0.4- 
0.2-- 



. 1 1 — I I I — — I 1 1 1 1 1 1 — I I 1 I— — I 1 — I — I 1 — I I I — I *■ 

S&P Nasd Cana Braz UK Irel Germ Autr Neth Swed Finl Spai Gree Indi SrLa Bang Japa HoKoTaiw SoKo Mala Indo Phil 

Figure 5: contributions of the stock market indices to eigenvectors e\, eis, and e23- Blue bars indicate positive values, 
and red bars correspond to negative values. 



One can see that the eigenvector corresponding to the largest eigenvalue is qualitatively different from the 
others. Nearly all markets (with the exception of Bangladesh and Indonesia) have positive representations. 
That is a compelling reason to believe that it represents a global market that is the result of the interactions 
of all local markets, or may also be the result of external news on the market as a whole. Figure 6 compares 
the time series of an index built using eigenvector 23 (in blue) with the world index calculated by the MS CI 
(Morgan Stanley Capital International), in red. Both indices are normalized so as to have mean two and 
standard deviation one. 



Indices 




Figure 6: time series of the market mode calculated using the eigenvector related with the largest eingenvalue of the 
correlation matrix (blue) plotted against the world index calculated by the MSCI (red). Both indices are normalized so 
as to have mean two and standard deviation one. 

In terms of portfolio theory, as stated by Markowitz' ideas |156j . |157j . the eigenvector corresponding to the 
largest eigenvalue represents the riskier portfolio one may build, as most of the indices vary in the same way. 
In constrast, some of the smaller eigenvectors represent portfolios with less risk, as, for example, eigenvector 
ei, which basically consists on "buying" S&P 500 (USA) and S&P TSX (Canada) and "short-selling" Nasdaq 
(USA), which are three very closely connected indices. Eigenvector e%§ corresponds to one of the eigenvectors 
that are within the region considered as noise, and should represent just a random combination of stock market 
indices. 

More differences between eigenvector e23 and the other eigenvectors can be seen if we build probability 
distribution of frequencies graphs for the twenty three eigenvalues. All distributions, except the one for eigen- 
vector e23, have average near zero and standard deviation around 0.21, while this is not the case for eigenvector 
e23- The elements of eigenvector e23 have mean 0.17 and standard deviation 0.13. 

Some recent works discussed how finite sized data and log-return distributions that are not Gaussian could 
affect the probability distribution of the eigenvalues of an empirical correlation matrix. Some of the results 
imply that the usual Marenko-Pastur distribution acquires a fat tail in the direction of the largest eigenvalue. 

A last analysis which shows the difference between the highest eigenvalues and the eigenvalues belonging 
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to the range associated with noise may be done using the so called Inverse Participation Ratio (IPR), 

N 



IPR k = £ (4) 4 , 



(7) 



where el is the i-th. element of eigenvector ej~, and N is the total number of eigenvectors. Its inverse gives the 
average number of stocks which contribute significantly to a portfolio built with such eigenvector. The next 
two figures show the IPR for the 23 eigenvectors, in ascending order from the left to the right (figure 7), and 
its inverse, PRk = 1/IPRk, for Participation Ratio (figure 8). 
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Figure 7: inverse participation ratio of the eigenvectors of the correlation matrix. 
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Figure 8: participation ratio of the eigenvectors of the correlation matrix. 

Note that, for eigenvector e23, the number of participating indices is larger than the average, which is about 
6. Most of the eigenvectors corresponding to noise fall around that average number, but this is not true for the 
eigenvectors corresponding to the lowest eigenvalues, which have a very small number of participating indices. 

One important result of this theoretical treatment is that the largest eigenvalue, associated with a market 
mode, is like another matrix that is added to the true correlation matrix of the log-returns. In order to study the 
remaining eigenvalues, one must first clean the empirical correlation matrix from the market mode. The process 
is known as single index model, and is widely used by theoreticians and practitioners of financial markets in 
order to remove the market mode of stocks negotiated in the same stock exchange |156| . This is done in order 
to study clusters of stocks that move together as blocks in stock markets. 

We now measure the average of the correlation matrices in a moving window of 30 days, changing one day at 
a time. The results are displayed in figure 9, where the average correlation is plotted together with the volatility 
of the market mode, which we consider here as the absolute value of St, where St is a linear combination of all 
indices with the elements of eigenvector e23 as the coefficients. The plot represents the average correlation of 
each window as a function of the last day of the window, so that events that occur after the date to which the 
average correlation is assigned do not influence its value. Volatility is in blue, and the average correlation is in 
red. 

Figure 10 shows the same information, but now both the average correlation and the volatility of the market 
mode are normalized so as to have mean two and standard deviation one. This is done in order to best compare 
both values, and it will be more useful for comparisons made for the other crises. 
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Figure 11 shows the average correlation and the average volatility of the market mode, both calculated in 
a running window of 30 days, and normalized so as to have mean two and standard deviation one. In this 
picture, the rise of valatility seems to be preceded by a rise in the correlation between international stock 
market indices, although that is not a conclusion that may be taken, since we are using averages here over a 
large period of time. 
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Figure 9: volatility of the market mode (blue) and average correlation (red) based on the log-returns for 1987, both 
calculated in a moving window of 30 days. 
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Figure 10: volatility of the market mode (blue) and average correlation (red) based on the log-returns for 1987, calculated 
in a moving window of 30 days and normalized so as to have mean two and standard deviation one. 
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Figure 11: average volatility of the market mode (blue) and average correlation (red) based on the log- returns for 1987, 
both calculated in a moving window of 30 days and normalized so as to have mean two and standard deviation one. 



It is quite clear that there is a strong correspondence between global market volatility and the correlation 
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of the market indices. The correlation between the two variables along this period is 0.62. One can also note 
that markets are much more correlated after the period of crisis, and this behavior tends to endure for some 
time after the crash [154], although one must take into account that the averaging procedure for the average 
correlation makes the curve smoother and thus decreasing less steeply. Figure 12 shows the evolution of the 
covariance between volatility and < C > in time, calculated in a moving window of 30 days, starting from 
02/12/1987 (the first day we assign an average correlation). A clear peak can be seen on the days of greatest 
volatility (we plot the covariance at the end of the time interval considered for each calculation). Although 
the covariance is influenced by the value of the volatility, so we expect to have large covariance when volatility 
is high, it has shown to be more efficient in determining periods of crisis than the correlation, that being the 
reason we are using it. 
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Figure 12: covariance between volatility and average correlation as a function of time. 



6 1998, Russian Crisis 

The Asian Financial Crisis, which ocurred in 1997, made the demand for raw materials fall worldwide, affecting 
Russia in particular, which is one of the major world exporters of commodities. With the war in Chechnya, and 
the transition to a capitalist economy, Russia showed signs of decline in its economy. By May, 1998, the fears 
concerning the Russian economy brought most of the world's financial markets down, since many countries had 
a good amount of money invested in that country. 

In order to analyze that crisis, we added to the previous indices the following: IPC from Mexico (Mexi), 
BCP Corp Costa Rica from Costa Rica (CoRi), Bermuda SX Index (Bermuda), Jamaica SX Market Index 
from Jamaica (Jama), MERVAL from Argentina (Arge), IPSA from Chile (Chil), IBVC from Venezuela (Vene), 
IGBVL from Peru (Peru), CAC 40 from France (Fran), SMI from Switzerland (Swit), FTSE MIB from Italy 
(Ital), BEL 20 from Belgium (Belg), OMX Copenhagen 20 from Denmark (Denm), OBX from Norway (Norw), 
OMX Iceland All-Share Index from Iceland (Icel), PSI 20 from Portugal (Port), PX from the Czech Republic 
(CzRe), PX from Slovakia (Slok), Budapest SX Index from Hungary (Hung), WIG from Poland (Pola), BET 10 
from Romania (Roma), OMXT from Estonia (Esto), PFTS from Ukraine (Ukra), MICEX from Russia (Russ), 
ISE National 100 from Turkey (Turk), TA 25 from Israel (Isra), BLOM from Lebanon (Leba), TASI from Saudi 
Arabia (SaAr), MSM 30 from Ohman (Ohma), Karachi 100 from Pakistan (Paki), SSE Composite from China 
(Chin), SET from Thailand (Thai), S&P/ASX 200 from Australia, CFG 25 from Morocco (Moro), EGX 30 
from Egypt (Egyp), Ghana All Share from Ghana (Ghan), NSE ASI from Nigeria (Nige), NSE 20 from Kenya 
(Keny), FTSE/JSE Africa All Shares from South Africa (SoAf), and SEMDEX from Mauritius (Maur). So, 
now we have a total of 63 indices, 5 from North America (if we include Bermuda), 2 from Central America and 
the Caribbean, 5 from South America, 24 from Europe, 2 from Eurasia, 17 from Asia, 1 from Oceania, and 

7 from Africa, where we are considering Russia and Turkey as part of Eurasia, for both countries are located 
in both continents. This offers a good degree of diversification, and includes Russia, which was of paramount 
importance in that particular crisis. 

Using the modified log-returns (JBJ based on the closing indices from 01/02/1998 to 12/30/1998, we built a 
63 x 63 correlation matrix between those. This matrix has average correlation < C >= 0.17, standard deviation 
a = 0.04, and is based on L = 257 days for the M = 63 indices, which gives Q = L/M = 257/63 « 4.079. 
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The upper and lower bounds of the eigenvalues of the Marenko-Pastur distribution ([3]) are 

A_ = 0.255 and A + = 2.235 . (8) 

The frequency distribution of the eigenvalues is displayed bellow (figure 13), plotted against the theoretical 
Marenko-Pastur distribution were it an infinite random matrix with mean zero and standard deviation 1. Figure 
14 shows the eigenvalues in order of magnitude, with the area corresponding to noise shaded. 

Note that the largest eigenvalue is completely out of scale. We also have several eigenvalues that are bellow 
the minimum theoretical eingenvalue and two other eigenvalues above the maximum theoretical eigenvalue. 

The next picture (figure 15) shows eigenvector e63, which corresponds to a combination of all indices in a 
market movement that explains about 36% of the collective movement of all indices. 
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Figure 13: frequency distribution of the eigenvalues of the Figure 14: eigenvalues in order of magnitude. The 
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Figure 15: contributions of the stock market indices to eigenvector e4s, corresponding to the largest eigenvalue of the 
correlation matrix. Blue bars indicate positive values, and red bars correspond to negative values. The indices are aligned 
in the following way: S&P, Nasd, Cana, Mexi, CoRi, Berm, Jama, Bra, Arg, Chil, Ven, Peru, UK, Irel, Fran, Germ, 
Swit, Autr, Ital, Belg, Neth, Swed, Denm, Finl, Norw, Icel, Spai, Port, Gree, CzRe, Slok, Hung, Pola, Roma, Esto, 
Ukra, Russ, Turk, Isra, Leba, SaAr, Ohma, Paki, Indi, SrLa, Bang, Japa, HoKo, Chin, Taiw, SoKo, Thai, Mala, Indo, 
Phil, Aust, Moro, Egyp, Ghan, Nige, Kcny, SoAf, Maur. 

Note that most indices have similar participations, with the USA and European indices appearing with the 
largest components for the eigenvector. The smallest participations, some of them with very small negative 
values, are the ones from Costa Rica, Bermuda, and Jamaica (Central America and the Caribbean), Iceland 
and Slovakia (Europe), all the Arab countries and most of the Southern Asia ones, China, and the African 
countries, with the exception of South Africa. 

Figure 16 shows the market volatility, together with the average correlation between the indices for 1998, 
using a running window of 70 days, and representing the average correlation of each window as a correlation of 
the last day of that window. The window has been enlarged due to the increase in the number of indices so as 
to avoid too much statistical noise. The volatility of the market mode is in blue, and the average correlation 
is in red. Both are normalized so as to have mean two and standard deviation one. This is done in order to 
better compare both measures. 

Note that the average correlation is high throughout the period, and it increases beginning in August, 1998, 
which is the start of the Russian crisis. The volatility of the market mode also grows higher during the same 
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period, although it presents some peaks prior to that time. As the market was unstable due to the Asian crisis 
of the previous year, that can be explained as well, although there was a drop in correlation between the world 
stock markets around April, 1998. 

Figure 17 shows the average volatility of the market mode (blue) and the average correlation, both nor- 
malized so as to have mean two and standard deviation one, and both calculated in a moving window of 70 
days. 

One can see that volatility and average volatility are correlated with the average correlation between the 
indices during the times of crisis. This dos not seem to be the case at the beginning of the year, when there 
was no crisis. 
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Figure 16: volatility of the market mode (blue) and average correlation (red) based on the log-returns for 1998, calculated 
in a moving window of 70 days and normalized so as to have mean two and standard deviation one. 
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Figure 17: average volatility of the market mode (blue) and average correlation (red) based on the log-returns for 1998, 
both calculated in a moving window of 70 days and normalized so as to have mean two and standard deviation one. 

The covariance between the volatility (not the average volatility) and the average covariance, in red, calcu- 
lated in a moving window of 30 days, is plotted in figure 18. One can verify that the covariance between them 
increases during the Russian crisis. 
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Figure 18: covariance between volatility (blue) and average correlation (red) as a function of time. 
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7 2001, Burst of the dot-com bubble and September 11 

On September, 11, 2001, the world was shocked, as the biggest terrorist atack in human history was perpe- 
trated against the USA. The death toll was close to 3,000, when two hijacked airplanes were thrown into the 
Twin Towers of the World Trade Center, in New York, one hit the Pentagon, in Virginia, and another fell in 
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Pennsylvania. Financial markets all over the world plummeted, in an uncertainty crisis that lasted a few days. 

On that same year, closer to March, there was the end of a financial bubble centered on internet-based 
companies, the so-called burst of the dot-com companies. That crash affected most markets in the world and 
is believed to be a result of an escalation of speculation with companies whose true values were much bellow 
the prices their stocks were being negotiated with. 

Here we analyze these two crises, one (September 11) which is a good exemple of a crisis which is caused by a 
completely exogenous cause, and the other (burst of the dot-com bubble) which is the result of high speculation 
on stock prices. For 2001, we use 53 indices, adding the KSE 100 from Pakistan (Paki), the Tunindex from 
Tunisia (Tuni), the SOFIX from Bulgaria (Bulg), the KASE from Kazakhstan (Kaza), and the NZSX 50 from 
New Zealand (NZ) to the ones already used for 1998. 

Using the modified log-returns © based on the closing indices from 01/02/2001 to 12/31/2001, we built 
a 79 x 79 correlation matrix between those. This matrix has an average correlation < C >= 0.11, standard 
deviation a = 0.03, and is based on L = 260 days for the M = 79 indices, which gives Q = L/M = 260/79 
3.29. 

The upper and lower bounds of the eigenvalues of the Marenko-Pastur distribution ([3]) are 

A_ = 0.295 and A+ = 2.122 . (9) 

The frequency distribution of the eigenvalues is displayed bellow (figure 19), plotted against the theoretical 
Marenko-Pastur distribution were it an infinite random matrix with mean zero and standard deviation 1. Figure 
20 shows the eigenvalues in order of magnitude. The region associated with noise is shaded. 
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Figure 19: frequency distribution of the eigenvalues of the Figure 20: eigenvalues in order of magnitude. The 
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The largest eigenvalue is once more completely out of scale. We also have several eigenvalues that are bellow 
the minimum theoretical eingenvalue. The next picture (figure 21) shows eigenvector 679, which corresponds 
to a combination of all indices in a market movement that explains about 19% of the collective movement of 
all indices. 
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Figure 21: contributions of the stock market indices to eigenvector 679, corresponding to the largest eigenvalue of the 
correlation matrix. Blue bars indicate positive values, and red bars correspond to negative values. The indices are aligned 
in the following way: S&P, Nasd, Cana, Mexi, Pana, CoRi, Berm, Jama, Braz, Arge, Chil, Vene, Peru, UK, Irel, Fran, 
Germ, Swit, Autr, Ita, Malt, Belg, Neth, Luxe, Swed, Denm, Finl, Norw, Icel, Spai, Port, Gree, CzRe, Slok, Hung, 
Pola, Roma, Bulg, Esto, Latv, Lith, Ukra, Russ, Kaza, Turk, Isra, Pale, Leba, Jord, SaAr, Qata, Ohma, Paki, Indi, 
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SrLa, Bang, Japa, HoKo, Chin, Mong, Taiw, SoKo, Thai, Viet, Mala, Sing, Indo, Phil, Aust, NeZe, Moro, Tuni, Egyp, 
Ghan, Nige, Keny, Bots, SoAf, Maur. 

Many of the indices have a very small participation, which amounts to no participation, due to possible error 
bars, and many others have almost no participation. The indices that have participation smaller than 0, 05 are 
the ones from Central America, Bermuda, Venezuela, Malta, Slovakia, Romania, Bulgaria, Latvia, Lithuania, 
Ukraine, Kazakhstan, all the Arab countries, with the exception of Saudi Arabia, Sri Lanka, Bangladesh, China, 
Mongolia, Vietnam, Malaysia, Indonesia, Philipines, and all the African countries, with the exception of South 
Africa. Given the size of those markets, this is within what was expected. The major contributions come 
from the North American countries, the major South American ones, most of Western and Central Europe, 
the Czech Republic, Hungary, Poland, Estonia, Russia, Israel, Hong Kong, South Korea, Singapore, and South 
Africa. 

Figure 22 shows the average correlation, calculated in a running window of 80 days, and the volatility of the 
market mode, both normalized so as to have mean 2 and standard deviation 1, since the correlation between 
both measures becomes more transparent in this framework. The normalized volatility is in blue, and the 
normalized average correlation is in red. 
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Figure 22: average volatility (blue) and average correlation (red) based on the log-returns for 2001, both calculated in a 
moving window of 80 days and normalized so as to have mean two and standard deviation one. 

The figure shows a great increase in volatility just after September 11, followed by an increase in average 
correlation between the world stock market indices. This is expected from a crisis that was completely exogenous 
to the financial markets. A similar increase of both volatility and average correlation occur close to the beginning 
of the year, related with the burst of the dot-com bubble. 

More illustrative is figure 23, which shows the average correlation and the average volatility, both calculated 
in a running window of 80 days, normalized so as to have mean 2 and standard deviation 1. The normalized 
average volatility is in blue, and the normalized average correlation is in red. 
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Figure 23: average volatility (blue) and average correlation (red) based on the log-returns for 2001, both calculated in a 
moving window of 80 days and normalized so as to have mean two and standard deviation one. 

The covariance between the volatility and the average correlation is plotted in figure 24, calculated in a 
moving window of 30 days. One can readily identify a peak around September 11, but no peak related with 
the burst of the dot-com bubble, which was not a precisely defined event in time. 
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Figure 24: covariance between volatility and average correlation as a function of time. 



8 2008, Subprime Mortgage Crisis 

The last large financial crisis began in 2007, reached its peak in 2008, and is happening until now. This crisis 
was triggered by the default of a large number of mortgages in the USA. Subprimes are loans to borrowers who 
have low credit scores. Most of them had a small initial interest rate, adjustable for future payments, which 
led to many home foreclosures after the rates climbed substantially. Meanwhile, the loans were transformed in 
pools that were then resold to interested investors. Since the returns of such investments were high, a financial 
bubble was created, inflating the subprime mortgage market until the defaults started to pop up. 

Because of their underestimation of risk, financial institutions worldwide lost trillions of dollars, and many 
of them declared banckrupcy. Because of that, credit lines tightened around the world, taking the financial 
crisis to the so called real economy. The world is yet to recover from this crisis, and many institutions are still 
to lose a good part of their assets in the following years. 

Here we analyze the year 2008, which is considered the time when the subprime crisis reached its peak, 
marked by events like the Lehman Brothers' announcement of banckrupcy, and the liquidation of three of the 
largest investment banks in the USA. In our research, we add now 13 indices to the ones we used for 2001: 
IGBC from Colombia (Colo), BELEX 15 from Serbia (Serb), CROBEX from Croatia (Croa), SBI TOP from 
Slovenia (Slov), SASE 10 from Bosnia and Herzegovina (BoHe), MOSTE from Montenegro (Mont), MBI 10 
from Macedonia (Mace), CSE from Cyprus (Cypr), Kwait SE Weigthed Index from Kwait (Kwai), Bahrain 
All Share Index from Bahrain (Baha), ADX General Index from the United Arab Emirates (UAE), DSEI from 
Tanzania (Tanz), and FTSE/Namibia Overall from Namibia (Nami). So, we use a total of 92 indices, 4 from 
North America, 2 from Central America, 2 from the islands of the Atlantic, 6 from South America, 35 from 
Europe, 2 from Eurasia, 28 from Asia, 2 from Oceania, and 11 from Africa. The sample became larger mainly 
because of the partition of the former Yugoslavia into many countries. 

Using the modified log-returns © based on the closing indices from 01/02/2008 to 12/31/2008, we built 
a 92 x 92 correlation matrix between those. This matrix has an average correlation < C >= 0.26, standard 
deviation a = 0.05, and is based on L = 253 days for the M = 92 indices, which gives Q = L/M = 256/92 ~ 
2.78. 

The upper and lower bounds of the eigenvalues of the Marenko-Pastur distribution ([3]) are 

A_ = 0.160 and A + = 2.558 . (10) 

The frequency distribution of the eigenvalues is displayed bellow (figure 25), plotted against the theoretical 
Marenko-Pastur distribution were it an infinite random matrix with mean zero and standard deviation 1. 
Figure 26 shows the eigenvalues of the correlation matrix in order of magnitude (the region associated with 
noise appears shaded). 



15 



1.2 



0.8 



0.4 



P(A) 




+ 



+ 



+ 



7A 



A 




+ 



+ 



+ 



2345678 31 
Figure 25: frequency distribution of the eigenvalues of the 
correlation matrix for 2008. The theoretical distribution 
is superimposed on it. 
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Note that the largest eigenvalue is even more out of scale than in previous crisis, what usually indicates a 
high level of correlation between the market indices and the presence of a powerful global market movement, 
although it is also influenced by the size of the sample of indices. We also have several eigenvalues that 
are bellow the minimum theoretical eingenvalue. The next picture (figure 27) shows eigenvector eg2, which 
corresponds to a combination of all indices in a market movement that explains about 34% of the collective 
movement of all indices. 
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Figure 27: contributions of the stock market indices to eigenvector eg2, corresponding to the largest eigenvalue of the 
correlation matrix. Blue bars indicate positive values, and red bars correspond to negative values. The indices are aligned 
in the following way: S&P, Nasd, Cana, Mexi, Pana, CoRi, Berm, Jama, Braz, Arge, Chil, Colo, Vene, Peru, UK, Irel, 
Fran, Germ, Swit, Autr, Ital, Malt, Belg, Neth, Luxe, Swed, Denm, Finl, Norw, Icel, Spai, Port, Gree, CzRe, Slok, 
Hung, Serb, Croa, Slov, BoHe, Mont, Mace, Pola, Roma, Bulg, Esto, Latv, Lith, Ukra, Russ, Kaza, Turk, Cypr, Isra, 
Pale, Leba, Jord, SaAr, Kwai, Bahr, Qata, UAE, Ohma, Paki, Indi, SrLa, Bang, Japa, HoKo, Chin, Mong, Taiw, 
SoKo, Thai, Viet, Mala, Sing, Indo, Phil, Aust, NeZe, Moro, Tuni, Egyp, Ghan, Nige, Keny, Tanz, Nami, Bots, SoAf, 
Maur. 



Indices with small negative contributions are those from Iceland, which suffered the effects of the crisis 
with greater impact than most of the other countries, Mongolia, Nigeria, Tanzania, and Botswana. Very 
small participations (less than 0.050) are related with the indices from Central America, the Atlantic Islands, 
Venezuela, Malta, Slovakia, Bosnia and Herzegovina, Montenegro, Kwait, Pakistan, Bangladesh, Vietnam, 
Ghana, and Kenya. Indices with strong participation (greater than 0.100) are those from Canada, Mexico, 
most South American countries, most of the European countries, Russia, Turkey, Cyprus, India, Japan, Hong 
Kong, Taiwan, South Korea, Thailand, Singapore, Indonesia, Philipines, Australia, Namibia, and South Africa. 
The surprise is the participations of the indices from the USA - S&P 500 (0.089), and Nasdaq (0.082), which 
are lower than expected. 

The following two figures show the relation between the average correlation and the volatility of the market 
index. Figure 28 shows the average correlation (in red) calculated in a moving window of 100 days, and the 
volatility of the market index (in blue), both normalized so as to have mean two and standard deviation one. 
Figure 29 shows the average correlation (in red) and the average volatility (in blue), both calculated in a moving 
window of 100 days, and normalized so as to have mean two and standard deviation one. 

One can see that the period of high volatility seems to be preceded by a period of high correlation between 
the stock markets of the world. Figure 30 shows the evolution of the covariance between the mean correlation 
and the mean volatility, calculated in a moving window of 30 days. 
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Figure 28: volatility (blue) and average correlation (red) based on the log- returns for 2008, calculated in a moving window 
of 100 days, both normalized so as to have mean two and standard deviation one. 
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Figure 29: average volatility (blue) and average correlation (red) based on the log-returns for 2008, both calculated in a 
moving window of 100 days and normalized so as to have mean two and standard deviation one. 
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Figure 30: covariance between volatility and average correlation as a function of time. 



day 



9 Normality tests for the correlation matrix 

In this section, we make tests in order to check wether the elements of the correlation matrix exhibit a normal 
or close to normal probability distribution or not. A first analysis of the data might lead us to believe it does. 
Observe the following graphics, with the skewness and kurtosis of the probability distribution obtained by 
considering the elements of correlations matrix (except those of the diagonal) calculated over moving windows 
(the size of the windows vary for each of the years that were considered) . 

For 1987 (figures 31 and 32), the size of the running window is 30 days. Note that, near the Black Monday, 
which occurred in October, kurtosis dropped substantially, what seems to imply that the distribution of the 
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coefficients of the correlation matrix approach that of a normal curve. 
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Figure 31: skewness of the correlation matrix during the year 1987 for correlation matrices calculated over a running 
window of 30 days. 
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Figure 32: kurtosis of the correlation matrix during the year 1987 for correlation matrices calculated over a running 
window of 30 days. 

For 1998, the running window has 70 days. The skewness (figure 33) remains nearly constant for most of 
the time, and the kurtosis (figure 34) of the same distribution stays near 3 during the same period. 
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Figure 33: skewness of the correlation matrix during the year 1998 for correlation matrices calculated over a running 
window of 70 days. 
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Figure 34: kurtosis of the correlation matrix during the year 1998 for correlation matrices calculated over a running 
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window of 70 days. 



The next two graphics show the skewness (figure 35) and the kurtosis (figure 36) for the elements of the 
correlation matrix, except its diagonal, calculated in a running window of 80 days for 2001. Both skewness and 
kurtosis present a peak in September 11, but otherwise remain nearly constant throughout the period. 
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Figure 35: skewness of the correlation matrix during the year 2001 for correlation matrices calculated over a running 
window of 80 days. 
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Figure 36: kurtosis of the correlation matrix during the year 2001 for correlation matrices calculated over a running 
window of 80 days. 

Figures 37 and 38 show the skewness and the kurtosis for the elements of the correlation matrix, except its 
diagonal, calculated in a running window of 100 days, for 2008. Note that the skewness becomes negative for 
the time after the beginning of the crisis, something that didn't happen in the previous cases. The kurtosis 
drops to values bellow 3 for the period of crisis. 
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Figure 37: skewness of the correlation matrix during the year 2008 for correlation matrices calculated over a running 
window of 100 days. 
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Figure 38: kurtosis of the correlation matrix during the year 2008 for correlation matrices calculated over a running 
window of 100 days. 



Since a perfect normal distribution would have skewness zero and kurtosis 3, we may see that the distribution 
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of the elements of the correlation matrix on international indices in periods of crisis are not normal, although in 
the case of 1987 and 2008, it seems to be the case. This assumption is contradicted if one plots the probability 
distributions of the correlation matrix every two months (figures 39, 40, 41, and 42). During the months of 
highest volatility of each crisis (October for 1987, August for 1998, September for 2001, and October for 2008), 
the probability distribution deviate somewhat from a normal distribution. 

One can see that the probability distributions for 2008 are less strongly peaked than for the other years, 
but this is mainly caused by the inclusion of a large number of weakly correlated indices. One can also notice 
that, in the months of crises, the average correlation grows, but the correlation gets more evenly distributed 
among the possible spectrum. 
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Figure 39: probability distributions of the correlation 
matrix calculated every two months in 1987. The 
probability distribution for the data for the whole 
year appears last, in red. 
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Figure 41: probability distributions of the correlation 
matrix calculated every two months in 2001. The 
probability distribution for the data for the whole 
year appears last, in red. 
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Figure 40: probability distributions of the correlation 
matrix calculated every two months in 1998. The 
probability distribution for the data for the whole 
year appears last, in red. 
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Figure 42: probability distributions of the correlation 
matrix calculated every two months in 2008. The 
probability distribution for the data for the whole 
year appears last, in red. 



Our claim that the probability distributions are far from Gaussian during periods of high volatility may be 
substantiated by using two tests for normality of those distributions. The Jarque-Bera test [158] is based on 
the formula 



JB 



N 



s 2 + 



(k - sy 



(ii) 



where N is the size of the sample, s is its skewness, and k is its kurtosis. The Lilliefors test [159], a variant of 
the Kolmogorov-Smirnov test, is based on the formula 



L = max\E(x) - N(x)\ 



(12) 



where E(x) is the cumulative distribution function estimated from the data and N(x) is the cumulative distri- 
bution function of a normal distribution with the same mean and standard deviation as the data. 

The Jarque-Bera test rejects the null hypothesis that the distribution is normal at the 5% significance level 
for all months of the years we have studied. The Lilliefors test rejects the null hypothesis that the distribution 
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is normal at the 5% significance level for all months except March/April and May/ June, 1987. When applied 
to the whole years of data, both tests strongly reject the hypothesis that the distribution of the correlation 
matrix is similar to a normal distribution. 

10 Gauging the results 

As we commented in the introduction of this article, one of the major concerns when dealing with data from stock 
markets all over the world is that most of them do not operate at the same hours. This leads to incorrections 
when one tries to study the correlations between them. Another source of concern is that sometimes the 
correlations between markets may not be measured correctly by the Pearson correlation coefficient, since it is 
better suited for linear correlation, which may not be the case. Other correlation coefficients, like Spearman's 
or Kendall's rank correlation coefficients, may capture relations which are not seen using Pearson's correlation. 

In order to gauge the effect of these two possible problems, we did two aditional analysis of the data. In 
the first one, we phased the data of Eastern markets (from Russia to the east) so that the data of Western 
stock markets were compared with data from the next day of Eastern markets. In the second one, we switched 
to Spearman's correlation whenever Pearson's correlation was used. We did all the calculations again for both 
cases and compared the results with the ones previously obtained. An account of the comparisons is given now 
for the four crises being considered. 

For 1987 with phased data, the average correlation becomes < C >= 0.15, slightly smaller than the value 

< C >= 0.16 for the unphased data. Using Spearman's correlation, we obtain < C >s= 0.07 (remember it 
is a different type of correlation, and so it should not be compared numerically with the Pearson correlation). 
For the phased data, the maximum eingenvalue, which was A max = 6.500, becomes A max = 6.135, and for 
Spearman's correlation, it becomes A max = 3.977. 

While for the original data Indonesia had a substantial negative participation in the eigenvector with the 
highest eigenvalue, no index has relevant negative participation for the phased data, and Taiwan and South 
Korea increase their participation, although Hong Kong decreases its own. For the eigenvector obtained with 
Spearman's correlation, Brazil, Finland, Bangladesh, and Taiwan acquire small negative participations, Sri 
Lanka and Indonesia maintaining their negative coefficients. 

For the phased data, there is nearly no change in the relations between average correlation and volatility, or 
between average correlation and average volatility, and the skewness and kurtosis of the probability distributions 
for the correlation matrix are also very similar. For the results obtained using Spearman's correlation, the 
agreement between the average correlation and average volatility is much greater for the data concerned with 
the begining of the crisis. 

For 1998, the average correlation < C >= 0.17 drops to < C >= 0.15 for the phased data and is given 
by < C >s= 0.16 for the data related with Spearman's correlation. The maximum eigenvalue goes from 
Amax = 16.897 to A max = 15.511 (phased data) and A max = 16.022 (Spearman's correlation). The participation 
of the Asian markets in the eigenvector corresponding to the largest eigenvalue grows for phased data and keeps 
essentially the same for Spearman's correlation. There are no substantial changes between average correlation 
and volatility and average volatility calculated in a moving window, nor in the skewness and kurtosis of 
the probability distribution of the off-diagonal elements of the correlation matrix, although for Spearman's 
correlation, the average correlation and the average volatility are slightly more connected. 

For 2001, the average correlation < C >= 0.11 remains < C >= 0.11 for the phased data, and it is 

< C >s= 0.07 for the data obtained using Spearman's correlation. The maximum eigenvalue goes from 
Amax = 15.284 to A max = 15.052 (phased data) and A max = 10.577 (Spearman's correlation). For the phased 
data, the number of participating Asian countries clearly grows, and the average participation, including those 
of some Western countries, also grows, but not substantially. For Spearman's correlation, the participations 
of indices in the market mode do not change substantially. There are no substantial changes to the skewness 
and kurtosis of the probability distribution of the off-diagonal elements of the correlation matrix if we use 
phased data or Spearman's correlation. For Spearman's correlation, the relation between average correlation 
and average volatility is even clearer. 

For 2008, the average correlation < C >= 0.26 drops to < C >= 0.21 for the phased data and is < C >s= 
0.22 for the data related with Spearman's correlation. The maximum eigenvalue goes from A max = 31.284 
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to A max = 26.761 (phased data) and A max = 26.587 (Spearman's correlation). The participation of Asian 
and African markets increase slightly in the eigenvector corresponding to the largest eigenvalue for the case 
of phased data. For Spearman's correlation, participations do not change significantly. The relation between 
average correlation and average volatility becomes stronger using phased data, and increases drastically for 
Spearman's correlation. There is nearly no change in the skewness and kurtosis of the probability distribution 
of the off-diagonal elements of the correlation matrix for phased data, but for Spearman's correlation the 
skewness and kurtosis curves become smoother. 

What we may conclude from this analysis is that the use of phased data gives occasional better results, but 
in general makes the average correlation between indices lower. So, we don't really have compelling reasons to 
use phased data. Now, for the calculations using Spearman's rank correlation, the agreement between average 
correlation and average volatility increases, sometimes drastically, as may be seen by comparing figures 43, 44, 
45, and 46 (shown next) with figures 11, 17, 23, and 29, respectively. 

< C > n ,< vol > n 
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Figure 43: average volatility (blue) and average correlation (red) based on the log-returns for 1987, calculated in a moving 
window of 30 days and normalized so as to have mean two and standard deviation one, using Spearman's rank correlation. 
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Figure 44: average volatility (blue) and average correlation (red) based on the log-returns for 1998, calculated in a moving 
window of 70 days and normalized so as to have mean two and standard deviation one, using Spearman's rank correlation. 
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Figure 45: average volatility (blue) and average correlation (red) based on the log- returns for 2001, calculated in a moving 
window of 80 days and normalized so as to have mean two and standard deviation one, using Spearman's rank correlation. 
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Figure 46: average volatility (blue) and average correlation (red) based on the log- returns for 2008, calculated in a moving 
window of 100 days and normalized so as to have mean two and standard deviation one, using Spearman's rank correlation. 

These four figures summarize what we have attempted here: to show that high correlation between world 
indices goes hand in hand with high volatility, possibly causing and definitely being caused by it. 
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11 Conclusion and future research 



Using the correlation matrices of the log-returns of a diversity of market indices during times of crisis, we showed 
that markets tend to behave similarly during times of high volatility. In the process, we verified the results 
obtained in a diversity of articles, but now applied to world financial market indices, and not to equities. Some 
of those results are that the probability distributions of the eigenvalues of the correlation matrices show peaks 
that are far off the maximum values predicted by Random Matrix Theory. Another result was the presence of 
certain combinations of indices that emulate a joint movement of most indices in what is called a market mode. 
An analysis of the probability distributions of the correlation matrices obtained show that those distributions 
are not normal and tend to flatten (low kurtosis) in times of crisis. 

We also showed that the relation of the average correlation and the average volatility (as calculated using the 
market mode) increases when one uses Spearman's rank correlation instead of Pearson's correlation, possibly 
highlighting nonlinear relations between them. The covariance between average correlation and the volatility 
of the market mode seems to be a good indicator of when periods of acute crises occur. 

Some direction for future research is to analyze how the techniques used in this work are modified if we 
consider that the frequency distributions of the log-returns are not Gaussian. Another topic that is being 
pursued is to study the hierarchies between the many indices and its evolution in times of crisis. For that, 
we shall use a distance measure based on the correlation between indices and build Minimum Spanning Trees 
and also Asset Trees in order to study cluster formation between indices |160| . Some of the results obtained 
here shall also be used in our studies of financial markets as coupled damped harmonic oscillators subject to 
stochastic perturbations |161j . 
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A Stock Market Indices 

The next table (table 2) shows the stock market indices we used, their original countries, the symbols we 
used for them in the main text, and their codes in Bloomberg. In the tables, we use "SX" as short for "Stock 
Exchange". Some of the indices changed names and/or method of calculation and are designated by the two 
names, prior to and after the changing date. 



Index 


Country 


Symbol 


Code in Bloomberg 


North America 




S&P 500 


United States of America 


S&P 


SPX 


Nasdaq Composite 


United States of America 


Nasd 


CCMP 


S&P/TSX Composite 


Canada 


Cana 


SPTSX 


IPC 


Mexico 


Mexi 


MEXBOL 


Central America 




Bolsa de Panama General 


Panama 


Pana 


BVPSBVPS 


BCT Corp Costa Rica 


Costa Rica 


CoRi 


CRSMBCT 


Caribbean 




Jamaica SX Market Index 


Jamaica 


Jama 


JMSMX 


British overseas territories 




Bermuda SX Index 


Bermuda 


Berm 


BSX 
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Index 


Country 


Symbol 


Code in Bloomberg 


South America 




Ibovespa 


Brazil 


Braz 


IBOV 


Merval 


Argentina 


Arge 


MERVAL 


IPSA 


Chile 


Chil 


IPSA 


IGBC 


Colombia 


Colo 


IGBC 


IBC 


Venezuela 


Vene 


IBVC 


IGBVL 


Peru 


Peru 


IGBVL 


Western and Central Europe 




FTSE 100 


I T n i r' o r\ \\ i n crr\ r\ tyi 
UIllLcQ xVlIl^QUIIl 


TTK 


u t\ /\ 


ISEQ 


J.1 tl till Li 


J.1 CI 


TSEO 


CAC 40 


r I clilCc 




PAP 


DAX 


KjVL IIlcLIiy 


Germ 


DAY 


SMI 


Switzerland 


own 


QA/TT 
olvll 


ATX 


A n stria 


Autr 


ATX 


FTSE MIB or MIB-30 


Italy 


Ital 


FTSEMIB 


Malta SX Index 


Malta 


Malt 


MALTEX 


BEL 20 


Belgium 


Belg 


BEL20 


AEX 


Netherlands 


Neth 


AEX 


Luxembourg LuxX 


Luxembourg 


Luxe 


LUXXX 


OMX Stockholm 30 


Sweden 


Swed 


OMX 


OMX Copenhagen 20 


Denmark 


Denm 


KFX 


OMX Helsinki 


Finland 


Finl 


HEX 


OBX 


Norway 


Norw 


OBX 


OMX Iceland All-Share Index 


Iceland 


Icel 


ICEXI 


IBEX 35 


Spain 


Spai 


IBEX 


PSI 20 


Portugal 


Port 


PSI20 


Athens SX General Index 


Greece 


Gree 


ASE 


Eastern Europe 




PX or PX50 


Czech Republic 


CzRe 


PX 


SAX 


Slovakia 


Slok 


SKSM 


Budapest SX Index 


Hungary 


Hung 


BUX 


BELEX 15 


Serbia 


Serb 


BELEX15 


CROBEX 


Croatia 


Croa 


CRO 


SBI TOP 


Slovenia 


Slov 


SBITOP 


SASE 10 


Bosnia and Herzegovina 


BoHe 


SASX10 


MOSTE 


Montenegro 


Mont 


MOSTE 


MBI 10 


Macedonia 


Mace 


MBI 


WIG 


Poland 


Pola 


WIG 


BET 


Romania 


Roma 


BET 


SOFIX 


Bulgaria 


Bulg 


SOFIX 


OMXT 


Estonia 


Esto 


TALSE 


OMXR 


Latvia 


Latv 


RIGSE 


OMXV 


Lithuania 


Lith 


VILSE 


PFTS 


Ukraine 


Ukra 


PFTS 


Eurasia 




MICEX 


Russia 


Russ 


INDEXCF 


ISE National 100 


Turkey 


Turk 


XU100 
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Index 


Country 


Symbol 


Code in Bloomberg 


Western and Central Asia 




KASE 






KZKAK 


CSE 


( , VTiVl 1 




CYSMMAPA 


Tel Aviv 25 




Isr£i 


TA-25 


Al Quds 


Palestine 


Pale 


PASISI 


BLOM 


T ,pKan on 




BLOM 

1 J 1_J V. J IV J. 


ASE General Index 


Jordan 


Jord 


JOSMGNFF 


TASI 


Saudi Arabia 


SaAr 


SASEIDX 


Kwait SE Weighted Index 


Kwait 


Kwai 


SECTMIND 


Bahrain All Share Index 


Bahrein 


Bahr 


BHSEASI 


QE or DSM 20 


Qatar 


Qata 


DSM 


ADX General Index 


United Arab Emirates 


UAE 


ADSMI 


MSM 30 


Ohman 


Ohma 


MSM30 


South Asia 




Karachi 100 


Pakistan 


Paki 


KSE100 


SENSEX 30 


India 


Indi 


SENSEX 


Colombo All-Share Index 


Sri Lanka 


SrLa 


CSEALL 


DSE General Index 


Bangladesh 


Bang 


DHAKA 


Asia-Pacific 




Nikkei 25 


Japan 


Japa 


NKY 


Hang Seng 


Hong Kong 


HoKo 


HSI 


Shangai SE Composite 


China 


Chin 


SHCOMP 


MSE TOP 20 


Mongolia 


Mong 


MSETOP 


TAIEX 


Taiwan 


Taiw 


TWSE 


KOSPI 


South Korea 


SoKo 


KOSPI 


SET 


Thailand 


Thai 


SET 


VN-Index 


Vietnam 


Viet 


VNINDEX 


KLCI 


Malaysia 


Mala 


FBMKLCI 


Straits Times 


Singapore 


Sing 


FSSTI 


Jakarta Composite Index 


Indonesia 


Indo 


JCI 


PSEi 


Philipines 


Phil 


PCOMP 


Oceania 




S&P/ASX 200 


Australia 


Aust 


AS51 


NZX 50 


New Zealand 


NeZe 


NZSE50FG 


Northern Africa 




CFG 25 


Morocco 


Moro 


MCSINDEX 


TUNINDEX 


Tunisia 


Tuni 


TUSISE 


EGX 30 


Egypt 


Egyp 


CASE 


Central and Southern Africa 




Ghana All Share Index 


Ghana 


Ghan 


GGSEGSE 


Nigeria SX All Share Index 


Nigeria 


Nige 


NGSEINDX 


NSE 20 


Kenya 


Keny 


KNSMIDX 


DSEI 


Tanzania 


Tanz 


DARSDSEI 


FTSE/Namibia Overall 


Namibia 


Nami 


FTN098 


Gaborone 


Botswana 


Bots 


BGSMDC 


FTSE/JSE Africa All Share 


South Africa 


SoAf 


JALSH 


SEMDEX 


Mauritius 


Maur 


SEMDEX 



Table 2: names, codes, and abreviations of the stock market indices used in this article. 
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